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ABSTRACT 

We present a detailed X-ray and optical study of a distant fossil system RX J1416.4+2315 
(z « 0.137), combining Chandra and XMM-Newton observations, optical photometry and 
spectroscopy. X-ray emitting hot gas imaged by both the Chandra and XMM-Newton shows a 
globally relaxed spatial distribution, supporting the idea that fossil groups are old galaxy sys- 
tems with no recent mergers. However, the diffuse X-ray emission shows signs of asymmetries 
in the core of the system. With a mean gas temperature of ~ 4.0 keV and total gravitational 
mass of 3. 1 X IQ^^Mq, within the virial radius, this is better described as a fossil galaxy clus- 
ter rather than a fossil group. The temperature profile shows no sign of a significant cooler 
core despite a cooling time dropping to 5 Gyr within the resolved core. We find a mass con- 
centration parameter C200 ~ 1 1 which is relatively high for a cluster of this mass, indicative 
of an early formation epoch. 

Using the spectroscopically identified cluster members we present the galaxy luminosity 
function for this fossil system. We measure the velocity dispersion of the galaxies to be ^ 
700 km s~^ based on 18 confirmed members. The dynamical mass is nearly twice the total 
gravitational mass derived from the X-ray analysis. The measured R-band mass-to-light ratio, 
within the virial radius, is ^ 440 M0/L0 which is not unusual for clusters of galaxies. The 
central giant elliptical galaxy has discy isophotes and spectral features typical of elliptical 
galaxies. 

Key words: galaxies: clusters: general - galaxies: elliptical - galaxies: haloes - intergalactic 
medium - X-ray: galaxies - X-rays: galaxies: clusters 



1 INTRODUCTION 

Though rich clusters of galaxies are the host for giant elliptical 
galaxies, their high velocity dispersion is less conducive to galaxy- 
galaxy interactions and mergers. In contrast, smaller galaxy sys- 
tems with low velocity dispersions are more merger efficient, and 
hence potential sites for the formation of bright and massive galax- 
ies. N-body simulations of clusters in a hierarchical cosmological 
model show that galaxy merging naturally produces a massive, cen- 
tral galaxy with surface brightness and velocity dispersion profiles 
similar to observed brightest cluster galaxies. The central galaxy 
forms through the merger of several massive galaxies along a fila- 
ment early in the cluster's history. 

If a merger efficient system such as a galaxy group remains 
undisturbed due to relative isolation from other massive systems, 
a giant elliptical galaxy will form as a result of the internal mul- 
tiple merger due to dynamical friction. The time-scale of such a 
process depends on the local density and is shorter at early for- 
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mation epochs. The resulting system should lack L* galaxies but 
the fainter end of the luminosity function should remain intact, as 
only a small fraction of the accreted mass is due to the dynamical 
friction of smaller galaxies over a Hubble time (Dubinski 1998j). 
With highly relaxed intra-group hot gas revealed by X-ray obser- 
vations and a gap of at least 2 magnitude between the first and sec- 
ond ranked galaxies jjones et al. Il2003h . galaxy groups known as 
"fossil groups" appear to be the observed products of the above 
mechanism ( Ponman et al. 1994). 

Ther e are about a do zen such system s iden- 
fified iPomnanetal. l 1 1994 IVikfalinin et al. I Il999t 
i Mulchaev & Zabludofj 19991: Ijones et al. I l200d : iRomer et 
l200CtlMatsushitjl200lLljones et al.ll2003h but few of them have 
been studie d in detail iKhosro^i^ii^ones & Pon man| |2004|: 
Sun et al. 2004). Based on the space density of fossils- ljonese^H 
(2003) estimate that fossil systems represent 8%-20% of all galaxy 
systems with the same X-ray luminosity, and are as numerous as 
poor and rich clusters combined. Therefore they account for a 
significant fraction of the galaxy systems in the universe, though 
observationally they are difficult to find, especially at higher 
redshifts. 
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About half a dozen fossils have been studied in the X- 
ray where they are most impressive. The nearest known fossil 
group NGC 6482, for instance, has shown interesting properties 
(Khosroshahi, Jones & Ponman 2004), such as an unusually high 
mass concentration and absence of a cooler core despite the very 
short cooling time. lSun et al. I^2004^ studied a more massive fossil 
using Chandra and XMM-Newton , and find a much smaller cool 
region than expected. If our arguments about the early formation 
epoch of the fossils are correct, these systems should be ideal lab- 
oratories to study cool cores. This is because of the absence of an 
obvious cooling flow removal mechanism, major cluster merging, 
and therefore other heating mechanisms can be looked into with 
less ambiguity. 

Fossils groups are equally important in their optical properties, 
as they provide a suitable environment to study the stellar properties 
of giant elliptical galaxies. Fossils are isolated and old systems en- 
abling us to make a direct comparison between the morphology of 
ellipticals grown in fossil groups and their counterparts at the core 
of galaxy clusters. This is further motivated by an earlier finding 
that there are enough fossils, from their space density, to provide 
the giant dominant ellipticals in clusters. There is a need for de- 
tailed optical studies to address various issues such as the mass-to- 
light ratio and luminosity function of galaxies, in addition to X-ray 
studies. The central galaxy warrants special attention, as the age of 
its stellar population can provide further clues to the cluster history, 
complementing the X-ray studies. 

This paper is organised as following: Section 2 briefly reviews 
the group properties and describes the observations, data reduction 
and preparations. The results from imaging and spectral X-ray anal- 
ysis are presented in section 3. Section 4 describes the distribution 
of mass. Results from optical and near-IR photometry, optical spec- 
troscopy and mass-to-light ratio are presented in section 5. A dis- 
cussion and concluding remarks are summarised in section 6. 

We adopt a cosmology with Ho = 70 km s~^ Mpc~^ and 
fim =0.3 with cosmological constant Qa = 0.7 throughout this 
paper. At the redshift of RX J 14 16-1-23 15 the luminosity and angu- 
lar diameter distances are 650 Mpc and 500 Mpc, respectively and 
1 arcsec = 2.44 kpc. 



2 OBSERVATION AND PREPARATION 
2.1 The group: 

The group RX J1416.4-I-2315 (hereafter J1416), z ^ 0.137, was 
first studied as part of a volume limited sample of spatially ex- 
tended X-ray sources compiled during the WARPS project (Wide 
Angle ROS AT Pointed Su rvey; Jscharfetal. 1997; Jones et al. 
Il99y ; Perlman et al. 2002)). Details of the initial sample selec- 
tion and the search leading to the volume lim ited sample of fos- 
sil groups can be found in I Jones et al. I i2003h . It is the most X- 
ray luminous source in the sample of 5 fossil groups studied by 
[jones et al. (2003), with a ROS AT estimated X-ray bolometric lu- 
minosity of 1.1 X lO^** erg s~^. In a deep optical image it would 
be classified as a galaxy group or poor cluster centred on an ex- 
tremely dominant, luminous giant elliptical galaxy. The very ex- 
tended X-ray emission is detected to a semi-major axis length of 
3.5 arcmin (512 /ifo^ kpc) in the ROS AT PSPC (Position Sensitive 
Proportional Counter) and is elongated in a direction similar to that 
of the central elliptical galaxy. 

T his galaxy cluster ap pears in the bright SHARC cluster sur- 
vey of iRomer et al. I i200Cl) . The redshift is confirmed, but it is not 
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Figure 1. Background spectra obtained from the source field (black) and 
blank sky dataset (grey). The blank sky spectrum was normalised by the 
flux of events detected outside the field of view. 



identified as a fossil system by iRomer et al. I i200(]l) . The target of 
this ROSAT field was a candidate cluster of galaxies, but at a high 
redshift (z> 0.3), and unrelated to the fossil system. 

The VLA FIRST catalog of radio sources gives a radio source 
at the centre of the giant elliptical galaxy with an integrated flux of 
3.39 mjy at 1.4 GHz. The radio source is extended on a scale of 4 
arcsec. 



2.2 Chandra observations 

The Chandra observation of this group was performed in Septem- 
ber 2001 (Obs ID 2024) for a total of 14.8 ks using Chandra ACIS- 
S. No major flares were detected therefore an exposure of 14.5 is 
used in the analysis. Though this observation helps to detect and 
position the point sources with a high spatial accuracy and gives 
enough counts to study the global properties such as the spatial dis- 
tribution and the core temperature of the hot gas, it does not provide 
enough counts for a more detailed analysis out to a large radius. 

We follow the same procedure as 

iKhosroshahi. Jones & Ponmanl (^2004) for the Chandra analy- 
sis except that we use a later version of the analysis tools (CIAO 
3.2 and CALDB 3.0). We therefore correct for the quantum 
efficiency degradation of ACIS, which increases with time. We 
use CTI corrected blank-sky background files in our imaging 
and spectral analysis as the ACIS-S3 chip is contaminated by the 
source emission due to the large extension of the X-ray halo. To 
study the soft diffuse X-ray emission, point sources were detected 
using CIAO "wavdetect" and removed. Their locations were then 
filled with surrounding background using a linear interpolation. A 
total of 20 point sources were detected. The brightest point source 
was an off-centre background AGN at 14:16:24.7, -(-23:16:14.9 
(J2000). 



2.3 XMM-Newton observations 

The large collecting area and field of view of XMM-Newton make 
the XMM-Newton observations attractive despite relatively poor 
spatial resolution compared to that of Chandra . The system was 
observed by XMM-Newton for 28ks in July 2003. The observation 
was split into two parts of 19ks and 9ks. The data were reduced and 
analysed with version 6.1 of the Science Analysis System (SAS) 
and with the most recent calibration database as of January 2005. 
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Figure 2. X-ray contours from the soft (0.3-2.0 keV) diffuse, point sources 
removed, X-ray emission detected by Chandra are overlaid on a 6x6 
arcmin^ size deep R-band optical image. 

The 19ks observation was rendered completely unusable by ex- 
tremely high background levels, though the second interval was 
mostly unaffected by high background periods. After cleaning, 
there were 8ks and 4ks of good time for the M PS and PN cameras 
respe ctively. Blank sky background datasets iRead and Ponmai] 
I2OO3I) were also prepared and cleaned in the same way. 

The reasonably small size of the source on the EPIC detec- 
tors meant that there was sufficient local background data to en- 
able useful comparisons to be made with the blank sky background 
fields to determine their suitability for data analysis purposes for 
this dataset. The spectrum of the background in the observation was 
compared to that in the blank sky data, after normalising them by 
the particle background flux measured from the events detected out- 
side the field of view of the telescopes. The resulting background 
spectra are plotted in Fig Q for the PN detector (the MOS spectra 
are similar). 

There is a significant excess in the source dataset at the soft 
end, due to higher Galactic X-ray background at the source posi- 
tion. The normalisation of the continuum above 2 keV is too high 
in the blank sky background, but the normalisation in the particle- 
induced fluorescent lines (1.5 and 8 keV) is correct. This means 
that the ratio of particle to hard X-ray background flux in the source 
data set is higher than that in the blank sky background. When the 
blank sky spectrum is normalised by the particle background flux, 
the other background components are overestimated, giving rise to 
a higher continuum level. 

This result suggests an alternative method of normalising the 
blank sky background; using the 2 — 7 keV flux. The upper limit is 
imposed to avoid another fluorescent line complex in the PN data. 
This method gives the correct > 2 keV continuum, but underesti- 
mates the fluorescent line flux. The effect of these two methods on 
temperature measurements is discussed later. 

2.4 Optical photometry and spectroscopy 

The imaging and spectroscopic observations of the group was 
performed using the observational facilities of the Issac-Newton 
Group of Telescopes (ING), Kitt-peak National Observatory 
(KPNO) and UK Infrared Telescope (UKIRT). The R-band image 
was obtained using the INT 2.5m wide field imager during a ser- 




Figure 3. X-ray soft (0.5-2 keV) diffuse emission contours from XMM- 
Newton overlaid on Chandra soft diffuse X-ray emission (left) and the 
R-band image of the central galaxy (right). Both the Chandra and XMM- 
Newton show asymmetries in the core. The white cross shows the location 
of the radio source from the FIRST catalog and matches the position of a 
point source detected at the centre of the cluster. 



vice time observation in August 2000. Unfortunately the conditions 
were not photometric, so further R-band imaging was obtained with 
the 8k mosaic camera at the University of Hawaii 2.2-m telescope 
in photometric conditions, and used to calibrate the original images. 

The spectroscopic observations were performed during a run 
to study several fossil groups using a multi-slit spectrograph on 
the KPN0-4m telescope on the 1 1-I3th March 2000. The Ritchey- 
Cretien spectrograph and KPC-lOA grating gave a dispersion of 
2.75 A/pixel over 3800-8500 A and with 1.8 arcsec slitlets a res- 
olution of 6A (FWHM) was achieved. Risley prisms compensated 
for atmospheric dispersion. Spectra were obtained through three 
slitmasks, with typically an hour exposure on each. The spectro- 
scopic data were reduced and analysed in the standard way using 
IRAF. 

The high quality near infrared images of the target were ob- 
tained in June 2004 using UIST/UKIRT Data reduction was per- 
formed using the Observatory Reduction and Acquisition Control 
(ORAC) data reduction tool and four 180 sec exposure mosaic im- 
ages were then combined to increase the signal-to-noise. 



3 X-RAY ANALYSIS AND RESULTS 
3.1 Spatial distribution 

Figure |2| shows the contours of soft (0.3-2.0 keV) diffuse X-ray 
emission from the Chandra observation on R-band optical image. 
The contours extended to 500 kpc along the semi-major axis of the 
X-ray emission. They appear relaxed and aligned with the stellar 
major axis of the giant elliptical galaxy. 

The XMM-Newton soft (0.5-2 keV) diffuse X-ray emission 
contours are also shown (Fig.|3j. The diffuse emission image was 
produced by adaptively smoothing an exposure-corrected image of 
the emission detected by the three XMM-Newton EPIC cameras to 
99% significance. Unlike Fig.|2| here the point sources are not ex- 
cluded. The off-centre AGN can be seen just North-West of the 
centre of the X-ray emission. The X-ray emission from the system 
at this energy range appears relaxed, overall, except in the core, 
where asymmetries are seen. The Chandra diffuse emission of the 
core shows a NE-SW elongation and a SW 'tail' (Fig.|2} which 
are also seen in the XMM-Newton contours (Fig|3). However, some 
of the other irregularities seen in the XMM-Newton image are not 
apparent in the Chandra data. 
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Figure 4. Radial surface brightness profile and the double /3-model fit to 
the soft energy (0.3-2.0 keV) X-ray emission from Chandra . The dotted 
curve is the /3-model fit to the inner region, ~100 kpc. The dark soUd curve 
is the /3-model to the outer region. The thick gray (green) curve represents 
the /3-model fit to the XMM-Newton surface brightness profile. 
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Figure 5. Temperature profile of J1416. The triangles (black) show the de- 
projected temperature profile from XMM-Newton analysis and the circles 
(red) is the same from Chandra analysis. 



3.1.1 Surface brightness distribution 

A radial surface brightness profile (Fig|4j was extracted from the 
inner 200 arcsec of the ACIS-S3 image in a soft energy range (0.3- 
2.0 keV) and fitted by a one-dimensional /3-model, 

E(r) = Eo[l + (-)r''+'^', (1) 

where ro and 13 are the core radius and index, and Eo the central 
surface brightness. This gives a core radius of ro= 28.2+443 arcsec 
and /3=0.47io 02 but is far from a good fit to the data. 

With the clear ellipticity in the two-dimensional X-ray dis- 
tribution and larger area covered by XMM-Newton EPIC we per- 
formed a two dimensional /3-model fit to the MOS 1 , MOS 2 and 
PN images at the same energy band. The core radius along the ma- 
jor axis is ro = 54 ± 4 arcsec, with /3 = 0.54 ± 0.02 and ellipticity 
of 0.34± 0.06. A two-dimensional /3-model fit to the Chandra data 
gives a ro = 36 ± 3 arcsec and /3 = 0.44 ± 0.01 with an ellipticity 
of 0.36. 

There is a large difference in the values of /3 and the core 
radius obtained from Chandra and XMM-Newton analysis. To 
find the root of the discrepancy and also because of a relatively 
poor fit to the Chandra data we performed a double-/3 fit by 
fitting inner and outer regions of the image each with a single 
/3-model. The parameters of these /3-models are, however, tied 
to give a continuous gas d ensity dKhosroshahi. Jones & PoimianI 
l2004l|Pratt & Amaudl2003h at the break radius, 100 kpc. We find 
roin= arcsec and /3i„=0.3lt^ |;^ and roo„t= 37.0t4 3 arc- 

sec and /3out=0.53lQ g2 with a reduced x^=l-32. Fig|4|shows this 
/3-model fit and the /3-model fit to the XMM-Newton surface bright- 
ness profile. 

3.2 Spatially resolved spectroscopy 

3.2.1 Chandra 

The ACIS spectra in three successive circular annuli were extracted 
and fitted with a de-projected (using XSPEC projct model) ab- 
sorbed APEC ( Smith et al. 2001) model. A fixed hydrogen column 
density of NH,gai = 2 x 10^" cm~^ was included in the model to 



account for Galactic absorption. The background was extracted for 
each annulus from the blank sky observations in the same detector 
regions as the source spectra. Each spectrum contains about 1000 
counts in bins of 20 counts each. Since we were not able to estimate 
the metallicity from the spectral fits, we fixed the metallicity to 0.5 
Zq, the mean value fiom XMM-Newton analysis on a similar scale. 

The unabsorbed flux from the central 200 arcsec is 1.22 x 
10^^^ erg s~^ cm~^, which corresponds to an X-ray luminosity, 
1.04x10** ergs/s (0.5-8.0 keV) and a bolometric luminosity of 
1.24x 10** ergs/s. The fitted column density was found to be con- 
sistent with the galactic level, though poorly constrained with an 
error of about 0.1 x 10^^ cm~^. 



3.2.2 XMM-Newton 

Spectra were extracted from 5 circular annular regions defined to 
contain > 1200 source counts. The effect of different background 
subtraction methods on the measured temperatures was investi- 
gated. A local background spectrum was extracted from an annular 
region free from source emission. A background spectrum was ex- 
tracted from the blank sky background in the same detector region 
as the source spectrum, and was corrected for the residuals between 
the local and blank sky backgrounds using the "double subtrac- 
tion" method described by (e.g. Arnaud et al. (2002)). These two 
methods were also repeated with the events weighted (using evig- 
weight) to correct for the telescope's vignetting. This weighting has 
the disadvantage of incorrectly boosting the non-vignetted particle 
background component. The two blank sky methods (ie, with and 
without weighting) were performed with the blank sky background 
normalised by either the flux of events detected outside the field of 
view of the telescopes or using the 2-7 keV flux. This gives a total 
of 6 background subtraction methods. 

The spectra from each detector were grouped to contain > 20 
counts per energy bin, and were fit simultaneously in the 0.4-7keV 
band by an absorbed APEC model. The absorbing column was 
fixed at the Galactic value as before, and the temperature, metal 
abundance and normalisation were free parameters in the fit. The 
temperature and metallicity of the PN and MOS spectral models 
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were tied together, but the normahsations were allowed to fit inde- 
pendently. 

The temperature profiles obtained with the non-weighted 
methods were found to be in excellent agreement, regardless 
of background subtraction method. However, the methods using 
weighted spectra were less consistent. As one might expect, the 
dispersion of temperatures increased with radius, as the back- 
ground becomes more important. The agreement between the non- 
weighted methods is excellent at all radii, while the temperatures 
measured with the (evig) weighted methods become inconsistent 
with each other, and the non-weighted methods at larger radii. This 
is perhaps unsurprising, as we have already seen that the non- 
vignetted particle background in the source dataset is relatively 
high fsection l23l . Based on these results, we chose to use the non- 
weighted, local background method for all of the XMM-Newton 
spectral analysis. 

The surface brightness distribution of J 14 16 is clearly ellip- 
tical, so the effect of this on the measured quantities such as the 
temperature profile is investigated by extracting spectra in ellipti- 
cal annuli. The ellipses were chosen to closely match the isophotes, 
and contain ~ 1200 source counts. The resulting temperature pro- 
file was consistent with that measured with circular annuli, and so 
the circular temperature profile is used throughout for consistency 
with our assumption of spherical symmetry in estimating the mass 
of the system. 

With the methodology thus established, the temperature pro- 
file was deprojected using the XSPEC projct model. Similarly we 
obtained a PSF-deconvolved temperature profile, however due to 
limited statistics we were unable to constrain the temperatures with 
the choice of bin sizes described above which is the largest to en- 
able us perform a mass profile analysis. We therefore did not use 
the PSF-deconvolved temperature profile for this study. 

3.3 Temperature and metallicity 

Fig |5| shows the deprojected temperature profile of the cluster 
based on various analysis methods out to 4 arcmin (~ 700 kpc). 
The XMM-Newton temperature agrees with the overall temperature 
found from Chandra observations of the central 130 kpc. Within 
the resolution of the data we find no significant drop in the temper- 
ature of the central annulus, in contrast to the observed temperature 
profil es of many relaxed cl usters of galaxies which present cool 
cores JVikhlininetal. l2005h . 

With limited counts from Chandra observations, our estima- 
tion of the metal abundances are based on the XMM-Newton ob- 
servations, however, the deprojected metallicities are very poorly 
constrained. The mean metallicity, within 400 kpc, is found to be 
0.23±0.076 solar metalhcity. 

3.4 Gas entropy and cooling time 

The entropy of the X-ray emitting gas is defined here as: 

S{r) = kT{r)/ne{rf^^keYcra, (2) 

where is the electron density. The shape of the entropy profile 
is of special interest since preheating models predict large isen- 
tropic cores. Using the /?-model fit to the surface brightness profile 
and the central gas density, obtained from the spectral fit, we can 
derive an analytical relation for the gas entropy providing that we 
model the temperature profile using an analytical function. Two en- 
tropy profiles shown in Fig|5|are obtained in this way by assuming 
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Figure 6. Entropy profile of the hot gas in J1416. The data points are ob- 
tained from emission measure for each annulus with its associated temper- 
ature. Solid and dashed lines ai'e the entropy profiles assuming polytropic 
and isothermal temperature profiles, respectively. 



isothermality (T « 4.0 keV) and a polytropic fit (7 = 1.1) to the 
temperature profile. In addition we obtain the gas entropy in each 
spherical shell (data points in Fig|6j by estimating the electron den- 
sity via the emission measure assuming a Raymond-Smith plasma 
emission model for each spherical shell and the L is the X-ray lumi- 
nosity from the shell with the volume V. The gas entropy at 0. lr-2oo 
(see section 3.5.1 for estimate of r2oo) is 5 *0.1 ~ 300 ± 50 keV 
cm^ similar to the average value obtained bv lPorunan et al. I i2003h 
for clusters with similar temperatures. No isentropic core is seen, 
though our spatial resolution is limited. This is consistent with other 
examples of group and cl uster scale s ystems re ported to show n o 
large isentropic cores (■Pratt & AmauA2003: Ponman et al. l2003h . 

The cooling time of the gas was calculated by dividing the 
total thermal energy of each shell by the associated emissivity, 
^coo! ^ I {ne+n^H)i<:T ^ whcrc £ is the bolomctric emissivity ob- 
tained by assuming a Raymond-Smith plasma emission model. The 
resulting profile is presented in Fig0 



3.5 Gas and total gravitational mass profiles 

The gas mass and the total gravitational mass profile can be derived 
from the gas density and the temperature profiles assuming that the 
gas is in hydrostatic equilibrium and is distributed with spherical 
symmetry. To support our assumptions we note that the distribution 
of the hot gas in this system is fairly symmetric and regular with the 
exception of a small region in the core. In addition, the comparison 
between the temperatures from circular and elliptical shells shows 
no significant difference in this case. The Chandra analysis of the 
X-ray surface brightness distribution shows that there is no signif- 
icant discrepancy in the results of the surface brightness profile by 
performing a one-dimensional analysis by selecting circular annuli 
or a two-dimensional fit. Finally Fabricant et al. ( 1984) show that 
the mass profile, from X-ray analysis, of an oblate or prolate system 
is similar to the one obtained under spherical symmetry assump- 
tion. 

The total gravitational mass is given by: 
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Figure 7. The estimated cooling time for each shell. The resolved core has 
a cooling time ~ 5 Gyr, yet we observe no cooler core in the temperature 
profile. 



Mgrav{< r) = - 



kT{r)r 



dlnp{r) dlnT{r) 



dlnr 



dlnr 



(3) 



where G and mp are the gravitational constant and proton mass and 

^l = 0.6. 

The temperature profile seen in Fig|5]is modelled by a poly- 
tropic model with (5 and core radius fixed to the values obtained 
from the surface brightness fit to the Chandra data giving 7 ~ 1.09. 
It is then straightforward to derive an analytical description for the 
total mass using the expressions for density under the polytropic 
assumption. However the error in the mass should be obtained 
indirectly . We use a Monte-Carlo metho d which is described in 
[Khosros hahi. Jones & Ponman (200 J) and lNeumann & Bohringeij 
>1995f) for our error estimation. In brief, the formula for the 
mass profile consists of two parts, the variation in gas density 
(din n/d In r) and in the temperature (r(r) and dlnT/dln r). For 
the latter we have generated 1000 physical temperature profiles 
with a temperature envelope defined by the observed errors in the 
temperature and linear interpolation. A physical temperature pro- 
file is one which guarantees a monotonically increasing mass with 
radius. In this simulation the density profile parameters were fixed 
at the values from the /3-model fit. In order to estimate the contri- 
bution of the density variation (d In n/d In r) to the mass error, the 
surface brightness profile at the location of each data point was fit- 
ted with d In n/d In r as the free parameter instead of j3. The error 
was estimated at the 68% confidence level using the error matrix 
provided by the fitting program. We checked the results by moni- 
toring the variation. The total error in the mass at each data point 
was then derived by adding the two uncertainties quadratically. 



3.5. 1 NFW mass profile and concentration 

Motivated by the results of numerous cosmological N-body sim- 
ulations and to enable us to make a direct comparison with mass 
profiles from other studies, we attempt to fit a NFW profile 
(Navarro et al. 1995) 
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Figure 8. Total gravitational mass and gas mass profiles of J1416. The upper 
solid (red) curve shows the gravitational mass profile. The dashed curve 
represents the NFW profile. The error envelope of the total gravitational 
mass is shown with vertical hashes. The lower solid (green) curve shows 
the gas mass profile. The concentration parameter is C200 = 11-2 it 4.5. 
The mass data points and errors are those estimated using MC simulations. 



to the total gravitational mass density. Integrated the mass profile 
for a spherical mass distribution leads to. 



Mtot{< r) = l&-KPm{rs)rl[rs ln(l + r/r^) - 



1 + r/r 



(5) 



where p,n{rs) is the density at rs. The mass concentration param- 
eter can then be defined as C200 = r2oo/ra. By definition, r2oo 
is the radius within which the mean gravitational mass density is 
200 times the critical density, pc{z). For our assumed cosmology 

(t7m=0.3,A=0.7), the critical density is pc(2) = |^ (1 + 2)^(1 + 

) ■ 

We fit the NFW mass profile (equation|5} to the five observed 
data points with errors obtained from MC simulation. The best fit 
profile has Vs = 97 kpc. r2oo is then calculated by extrapolating 
the profile to 200pc. This gives r2oo = 1.22 ± 0.06 Mpc. The 
extrapolated total mass, using the best fit NFW profile, at r2oo is 
M200 = 3.1(±1.0) X 10^* M0. 

To obtain the concentration of the dark matter we take into 
account the contribution of the stellar mass of the central dom- 
inant galaxy, for which we have the R-band luminosity and as- 
sume a stellar M/Lr ~ 5 Mq/Lq from the recent fundamental 
plane studies ( CaoDellariet al. 2005). This gives — 108 kpc and 
therefore a concentration parameter C200 = 11-2 ± 4.5, which is 
high compa red to the values of C200 for poor clusters studied by 
IPratt & Am aud (2005). Implications of this high mass concentra- 
tion is discussed in section 5. 

The integrated gas mass can be derived by integrating the gas 
density. Fig|8|shows the gas mass profile. The total gas mass within 
r200 is ~ 4 X 10^^ M0. With the extrapolated total gravitational 
mass within the same radius ~ 3.1 x 10^'' M0, the gas fraction is 
calculated to be ~ 0.13. 
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Figure 9. Galaxy luminosity function of J1416. The dark data points and 
the associated en'or bai's show the luminosity function of the galaxies using 
a statistical subtraction of the background using two regions outside the 
virial radius. The curve (blue) is the best fit Schecter function to these data 
points. The open circles (red) represent the luminosity function of the the 
confirmed cluster members. 

4 OPTICAL PHOTOMETRY AND SPECTROSCOPY 

To estimate the completeness of our spectroscopic observation we 
extract galaxies down to i? = 19.5 mag using Sextractor. This 
gives a total of 48 sources wi thin 0.5r2oo (the r adius used in the 
definition of fossil systems bv ljones et al. I i2003h ). of which 9 are 
classified as stars. From the total of 39 galaxies, only 27 were ob- 
served, of which 18 were cluster members with a mean redshift of 
0.137 and calculated galaxy velocity dispersion of cr = 700 ± 120 
km s^^. Some of the galaxies (8 in total) have been missed due 
to the slit configuration. Four galaxies were know to be fore- 
ground galaxies prior to this observation. We estimate the com- 
pleteness of the sample to be 77 percent within 0.5r2oo and down 
to niR = 19.5. Based on this calibration, the magnitude of the cen- 
tral galaxy is 15.17 and the second brightest galaxy in the cluster is 
m_R = 17.7. None of the galaxies missed due to the slit configu- 
ration were brighter th an R = 17.7, thus t he system does meet the 
formal fossil criteria o f IJones et al. N2003h . 



siderably shallower. There is an inconsistency in the slopes of the 
galaxy luminosity functions of the confirmed members and that ob- 
tained from statistical background subtraction, which could be due 
to the fact that the statistical background subtraction is performed 
based on the galaxy classification, using SExtractor, in the absence 
of galaxy colours. 

For RX J1552.2+2013, another fossil with a simila r mass 
to J1416, Mendes de Oliveira. Cvpriano & Sodre Jr.l ^OO?) find a 
slope of a = — 0.77±0.37 for the spectroscopically selected galax- 
ies and a = — 0.64±0.30 for the photomet rically selected gala xies. 
From the study of the GEMS X-ray groups- lMOes et al. liioO? ) find 
a = —1.0. They also report a dip in the luminosity function which 
is more prominent in poor groups and is argued to be the result of 
galaxy mergers in groups. A similar mechanism is thought to be 
the origin of the large gap in fossils' luminosity function. This is 
discussed in section 5. 

4.2 Mass to light ratio 

Having the total gravitational mass and the optical luminosity of 
the cluster we are in a position to estimate the mass-to-light ratio 
of the system. 

The total R-band luminosity of the cluster is measured to be 
jqII.85 accumulating the luminosity of the photometrically 

identified cluster members within r2oo- Thus the total mass-to-light 
ratio within 7-200 is ~ 440 Mq/Lq. Converting this to the mass- 
to-light ratio in B-band we find M/Lb ~ 580 Mq/Lq. This 
M / L, derived assuming B-R=1.5, is an upper limit because all the 
cluster galaxies are assumed to be early-types. Under a similar as- 
sumption Mendes de Oliveira, Cvoriano & Sodre Jr. (2005) find a 
mass-to-light ratio of M/Lb ~ 750 A'Iq/Lq for another mas- 
sive fossil based on a dy namical mass estimate. These mass to light 
ratios, ours and that of iMendes de Oliveira. Cvpriano & Sodre JrJ 
i2005h . are relatively high compared to the results of dis- 
tant clusters M/L fl ~ 2IOM0/L0 ( Carlberg et al. 199^ 
lGirardietal.ll2002h . However, higher v alues of mass to light ra- 
tio are al so reported jMohr et al. Ill993). In c omparison to other 
fossils I Vikhlin in et al. Ill999l: lYoshioka et al. ||20 04) these fossils 
and NGC 6482 /Kh osroshahi. Jones & Pon m£m2Q04) have lower 
mass-to-light ratio and hence difficult to argue whether or not fos- 
sils have, in general, higher values of M/L. 



4.1 Luminosity function 

The luminosity function (LF) of galaxies in this fossil cluster is 
obtained by selecting the galaxies down to AIr ~-18.5. The cu- 
mulative LF is shown in Fig |9| after statistical subtraction of the 
background galaxies extracted from two fields around the source 
just outside the r2oo of the cluster. In this figure we show the lumi- 
nosity function of the spectroscopically confirmed cluster members 
(open red circles) and the luminosity function of galaxies (filled 
circles) after a statistical subtraction of the field galaxies and the 
associated error bars. 

Fitting a Schecter function to the LF of confirmed 
members, we obtain a = —1.23 ± 0.28 and — 
—20.40 ± 0.22. An independent study of a higher quality 
galaxy sample for this fossil cluster res ults in a similar slope 
iCvpriano. Mendes de Oliveira & Sodre Jr . 2006). We obtain a — 
-0.61 ± 0.20 and = -20.55 ± 0.40 when we fit a Schecter 
function to the LF of photometrically selected galaxies (solid curve 
in the Fig|9j. As this figure shows, this luminosity function is con- 



4.3 The central galaxy 

To study the stellar light distribution of the central giant elliptical 
galaxy we performed an isophotal analysis using the high resolu- 
tion UIST K-band images of the system. The isophotal analysis us- 
ing the IRAF task 'ellipse' shows a discy nature for the stellar light 
at an outer radius of r > 20 kpc as shown in figure Fig[T0| The 
ellipticity of the galaxy increases to a value of 0.43, overtaking the 
ellipt icity of the X-ray distribution as expected iBuote & Canizare3 
ri99w). To find out if this is a pure elliptical galaxy we fitted a 
Sersic profile (r^''") to the radial surface brightness profile us- 
ing a two dimensional decomposition technique which allows for 
a possible disc component, with an exponential profile. We find 
no disc contribution to the total galaxy light. We find a Sersic in- 
dex of n ~ 3.0 and a half light radius of ~ 20 kpc. The giant 
elliptical galaxies in galaxy clusters often show boxy isophotes 
(Bender. Burstein & Faber 1992) while having larger values of n 
i Khosroshahi et al. ..2004.) . Further discussion should be based on 
a well defined sample of dominant galaxies in fossils, which is the 
subject of a separate study. 
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Figure 10. Isophotal analysis of the central galaxy in K-band. Surface 
brightness radial profile (top-left), ellipticity profile (top-right), Cosine co- 
efficient B4 (bottom-left) and the stellar spectrum of the giant elliptical 
galaxy are presented. 

The spectrum of the central galaxy is similar to that of a typical 
elliptical galaxy, with H(5 and Mg lines marked in Fig llOl 



5 DISCUSSION AND CONCLUSIONS 

Using the high angular resolution of the Chandra and higher sen- 
sitivity of the XMM-Newton , combined with optical imaging and 
spectroscopy and near IR imaging, we studied the most massive 
and hot fossil galaxy system, which also hosts the the most lumi- 
nous giant elliptical galaxy observed in fossils. 

We confirm the result of an earlier ROSAT study that RX 
J1 416.4-I-2315 meets t he requirements of being a "fossil" discussed 
in Ijones et al. I i2003h . These are the minimum X-ray luminosity 
of lO"*^ erg s-i and a minimum of two magnitude difference in 
R-band luminosity between the two brightest galaxy belonging to 
the system (within O.Srvir)- The examination was necessary as the 
gas temperature in the earlier study was underestimated by ~ 40% 
resulting in a smaller virial radius within half of which we had to 
examine the above optical criteria. 

The global X-ray distribution is relaxed, except a small re- 
gion at the core, indicating an absence of any recent group- 
group merger and major galaxy-galaxy merger indicative of an 
early formation epoch. The central dominant galaxy is a gi- 
ant elliptical galaxy aligned with the X-ray emission and hence 
the dark matter distribution. The alignment which i s also seen 
in NGC 6482 jKhosroshahi. Jones & PonmarJ l2004b and ESO 
3060170 jSunetal.lbOoir is an interesting feature which is 
known to be the consequence of an anisotropic collapse, such as 
along a filament. The alignment effect only works for the cen- 
tral g iant elliptical g alaxy and is seen in the N-body simula- 
tions jDubinskilll99a) . Recent hydrodynamical cosmological sim- 
ulations, to investigate the effect of baryonic dissipation on halo 
shapes I Kazantzidis, Zentner & Nagai 2005), show that there is 
also a strong correlation between dark matter halo shape and stellar 
remnant morphology in collisionless mergers. 

Despite the high degree of relaxation and symmetry, there are 
signs of asymmetry at the core of the system. These are seen in 



the X-ray soft emission within the core. This could be a result of 
AGN activity. VLA FIRST catalog gives the integrated flux of the 
central AGN in the giant elliptical galaxy as 3.39 mjy at 1.4 GHz, 
with a spatial size of 4 arcsec. Radio observations of the system are 
currently being performed to investigate the detailed morphology 
of the radio emission and the status of the AGN duty cycle. 

If fossil groups are old systems, as argued, they should provide 
ideal environments for the formation of cool cores, due to the ab- 
sence of any recent major merger which can remove the cool core. 
The temperature profile shows no significantly cooler core, sug- 
gesting the existence of a heating mechanism. Given the limited 
statistics and the XMM-Newton point spread function, the presence 
of a small (< 50 kpc) cooler region cannot be ruled out. However 
the absence of cool core is also found in the case of the neares t 
observed fossil, NGC 6482 jKhosroshahi. Jones & Ponmanll2004 
where the core is resolved with 1 kpc resolution. In addition ESO 
3060170 iSunetal. 2004) shows a cool core much smaller than 
the expected cooling region for a system of few 10^^ A/q . As men- 
tioned above there is evidence of AGN activity in the core of RX 
J1416-I-2315, similarly to ESO 3060170. 

This study shows that the X-ray temperature from a previ- 
ous ROSAT study was underestimated by about half its value. It 
is therefore important to study the effects of these changes on the 
scaling relations involving X-ray quantities. This is the subject of a 
separate study concentrating on the X-ray scaling relations of fossil 
galaxy groups. 

The initial argument for the early formation of fossils re- 
lied on the relaxed X-ray morphology, absence of recent merg- 
ers and the large gap (2 mag) in the optical luminosity function 
of their constituent galaxies. Their mass distribution, and in par- 
ticular their mass concentration parameter, can pr ovide more solid 
and independent evidence for their early formation. lWechsler et al. I 
1 2002) show that halos which have not had a major merger since 
z = 2 are more concentrated in comparison to systems which 
had experienced a major merger since then. We measure C200 
and find that the mass concentration for t his system is relativel y 
high compared to clusters of similar mass jPratt & ArnauSbOOSh . 
This is excluding the stellar mass contribution to the total gravi- 
tational mass from the central giant elliptical galaxy. We have al- 
ready repo rted a very high mass co ncentration paramete r for NGC 
6482 I Kho sroshahi. Jones & Ponm an 2004). In addition. lSun et al.1 
1 2004J find a high mass concentration (~ 8.9) when they fit a NEW 
profile for ESO 3060170, another massive fossil. If fossils are old 
systems, their core must have formed in an early stages when the 
density of matter was higher, resulting in formation of a dense and 
compact core whilst later accretion of mass builds an extended halo. 
Relatively higher mass concentration of low mass systems com- 
pared to more massiv e halos are seen in cosmological simulations 
teuUocketal. 1200 ih . 

The central galaxy in this cluster has discy isophotes despite 
the fact that majority of giant elliptical galaxies at the cores of clus - 
ters and groups show boxy isophotes teor mendv & BendeJl99^ . 
There are other examples of discy isophotes in fossils' central 
galaxies. NGC 6482 also shows highly discy isophotes. The orig- 
inal fossil, RX J1340-f40I7 has pure elliptical isophotes. We are 
currently studying the possibility that this, non-boxy isophotes, is a 
generic property of giant ellipticals at the core of fossil systems. 

A more accurate galaxy luminosity function, compared to 
the first published for a fossil (Jones etal 2001), is presented 
in this study. The confirmed members form a LF which is 
steeper than most of the galaxy luminosity functions presented 
for the groups. We find a = -1.23 ± 0.28 and A/^ = 
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-20.40 ± 0.22, by fitting the Schecter function to the LF of 
the spectroscopically confirmed member galaxies. This agrees 
well with the results of an independent study of this cluster by 
[Cvp riano. Mendes de O liveira & Sodre Jr. ( 2006). 

iMiles et al 1120041) found, in their study of GEMS groups, 
an interesting feature referred to as "dip" in the optical luminos- 
ity function of dim (low velocity dispersion) X-ray groups com- 
pared to bright (high velocity dispersion) X-ray groups interpreted 
as being produced as a result of more efficient galaxy merger in 
low velocity dispersion groups. Similarly, the main argument be- 
hind the lack of bright galaxies in fossils is the merger of L* galax- 
ies, but fossils are in gen eral more X-ray lu minous than the dim 
X-ray groups studied by ( Mil es et al. Il2004h . The velocity disper- 
sion of galaxies in this fossil cluster is cr « 700 km s~^ which 
is con sistent with the value ex pected from the observed cluster 
Lx-a iMahdavi & Geller'200l|) given the X-ray luminosity of this 
cluster. Sim ilarly J1552.2+2 013, lies on the sam e Lx-<y relation 
iMendes de Oliveira. Cvpriano & Sodre Jrl2005h . though the scat- 
ter in the Lx-<J relation is large. 

An especially interesting feature of RX J1416.4+2315 is the 
large size and mass of the system. Dynamical friction, the key pro- 
cess believed to be responsible for the orbit decay which allows the 
major galaxies to merge within fossil groups, is much less effective 
at high velocities, raising the question of whether the existence of 
such a fossil cluster poses a problem for th e whole model of fo ssil 
system formation through galaxy merging. Ijones et alH i2000h es- 
timated the time-scale of dynamical friction for the original fossil 
group, J1340.6-I-4018, and showed that L* galaxies, initially in cir- 
cular orbits at half the virial radius and circular velocity of \/2a, 
should merge into the central galaxy in 4.5 x lO' yr. The subject 
of this study, J1416, is more massive and has a radial velocity dis- 
persion, (7 almost twice that of J1340.6+4018. It also has a larger 
virial radius. The time scale for orbital decay of a body of mass M 
by dynamical friction is proportional to ^-j^, teinnev & Tremain3 
[l987), which means that a L* galaxy in this system will require a 
dynamical friction time-scale larger than a Hubble time to merge. 
We note, however, that this estimate is based on initial conditions 
in which the galaxy is part of a virialised system and rotates on 
a circular orbit around the centre. If instead, galaxies fall into the 
developing cluster along a filament, the loss of the angular momen- 
tum could occur in a much shorter time, due to the small pericen- 
tre radius of their orbits. The X-ray emission, and hence the dark 
matter distribution in J1416, is highly elongated, implying a very 
anisotropic velocity dispersion tensor, and supporting the idea of 
collapse along a dominant filament. 

Recent simulations of the formation of fossil systems also sug- 
gest th at they can be formed with high masses. IP'On ghia et al I 
yOOjj) find that fossil groups have already assembled half of their 
final mass at z~l, and subsequently they typically grow by minor 
merging only. The early assembly of fossils groups leaves sufficient 
time for L* galaxies to merge into the central galaxy by dynamical 
friction, resulting in the large magnitude gap. 
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